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l/" - ) ■ Abstract 

i | A grand canonical ensemble of interacting bosons is considered. The zero 

temperature phase diagram is evaluated from the mean-field approximation 

^ i ■ of the functional integral. Three phases are found: a superfluid, a normal 

<^ , fluid and a Mott insulator. 

00 ■ 

Quantum systems with many-body interactions are an important application of 
Q\ . functional integrals. As an example a model of interacting bosons is considered. It 

+^ \ is denned by the Hamiltonian 

£ v 

H=- Y, JryH$r> - uYH®r + -Y®l®M®r, (1) 

£H ■ <r,r'> r r 

o ' 

where $J (<& r ) creates (annihilates) a boson at the lattice site r. The first term 
describes bosons hopping from site r' to a nearest neighbor site r. u is a single 
^ . particle potential and V the coupling constant of an on-site interaction. 

Following the standard procedure (e.g. Ref. 0) the corresponding functional 
integral for a grand canonical ensemble of bosons reads 



exp(-S)D[0,0*] (2) 



with the action 



S = A E 7T^ r ' *) *) ~ ^ " A )] + A E l H lP(r, t), 0(r, t - A)]. (3) 

We have introduced a discrete imaginary time with interval length A and a complex 
field with periodic boundaries in time <p(t = 0, r) = 0(t = f3H,r). The parameter 
plays in (|2|) the role of the chemical potential. 
The local part of the Hamiltonian in S can be approximated (f>*(r, t)<p(r,t — A) w 
<j)(r,t)(j)*(r,t) = \<f)(r, t)\ 2 , neglecting terms of 0(A 2 ) in S.JlJ The interaction term 
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has a positive coupling constant V (repulsion). Thermal fluctuations are switched 
off by taking implicitly the limit (3 — > oo. A random path expansion ||] can be 
created by extracting the local part from the action 

5 = -A^[(/i-l/AK-^.], (4) 

where n x = \<f) x \ 2 and H = 1. The expansion in terms of S\ = S — Sq generates 
locally powers of the complex field <p x and <f>*. Since Sq does not depend on the 
phase of <p x we obtain from the integration w.r.t. (f) nonzero quantities only if <j) x 
and 4>* appear with the same power at each point x = (t, r) 

J ftp* exp {A[(/i - 1/A)n x - = gi J L>l *. (5) 

Thus, the random path expansion of the partition function leads to a grand-canonical 
ensemble of Bose world lines (BWLs), going from t — to t — [3, which are directed 
along the t-axis. The interaction of the bosons corresponds to the coefficient g\ at x 
which counts I visiting bosons at r at the time t. The statistics of the BWLs simplifies 
essentially if we apply a hard-core approximation by introducing the restriction that 
a point x can be visited at most by one boson 

gi = for I > 2. (6) 

This condition projects onto boson realizations with hard-core interaction at each 
point. The contributions to the expansion are empty points with weight go = ( and 
world line elements which connect points x and x' with statistical weight 

!Ag 1 J r y = a for x' = x + e'^ 
1 for x' = x + Ae t ■ (7) 

otherwise 

are the unit vectors in the (i-dimensional space, e t is the unit vector in time and 
e' = e M + Ae t . This model has only two independent parameters, a and (. We may 
write for the partition function of a system on a lattice with M sites 

£ P ({rM}n), (8) 

where P({rj(t)} n ) is the weight of a given configuration {rj(t)} n of n BWLs 

P({r 3 ( t )}n)= II W *J- ( 9 ) 

x,x'E{rj(t)} n 

I n are all configurations of BWLs going from t = to t = /3 which are allowed by 
the hard-core condition. It has been shown that the partition function of hard-core 
bosons (H) can be rewritten for dimensions d > 1 in terms of interacting fermions 

^ferm = ^{w x , x > + C^x'^l^l'^l' , (10) 



2 



with Grassmann fields t/>£. The hard-core interaction is now represented by the 
exclusion (Pauli) principle of the fermions. Especially in d = 1, the hard-core 
bosons cannot exchange their positions. This implies that an additional interaction 
between the fermions is not necessary, and the hard-core bosons are equivalent to 
non-interacting fermions 0. 

The hard-core boson model has been generalized to a model of bosons with iV 
colors 0] or a model with N levels at each point [[J . In the latter the bosons can oc- 
cupy the N levels statistically. In both cases the limit N — > oo can be solved exactly. 
The formal reason is a formation of composite fields from pairs of fermions. For the 
iV color bosons these are nonlocal fields and in the case of the N bosonic levels these 
are local fields. The composite fields do not depend on the number N but iV appears 
as a prefactor of the action. Therefore, a saddle point integration can be performed 
for N — > oo in both cases. The results of the calculations of Refs. @, [|] will be used 
in the following to derive the phase diagram of hard-core bosons at zero temperature. 

a) local order parameter. Superfluid (SF) 

The effective functional integral of the N level model of Ref. || is related to a 
local complex order parameter field ((p x ,Xx)- The corresponding partition function 
reads 

Z loc = f exp { - N[(<p, V) + (x, X*) ~ E lQ g[C + (f* + {ft + }v[f, x], 

J X 

(11) 

where V\ is the matrix w + 1. There are two mean-field solutions for the saddle 
point equation Ss\ oc = 0: a trivial one with (p — x — and a nontrivial one with 
if = —ix = — (. The stability of the two solutions is found from the fluctua- 
tions around the saddle point ||: the trivial solution is stable for ( > a 4 and the 
nontrivial one for < ( < a 4 . Thus the critical point is ( c = a 4 . The mean-field 
approximation of the action is — log £ for the trivial solution and a 2 — 4 log a — (/a 2 
for the nontrivial solution. The fluctuations are gapless for < ( < a 4 due to Gold- 
stone bosons created by a spontaneously broken £7(1) symmetry. The corresponding 
phase is a superfluid. This result can also be obtained if a Ginzburg-Landau ex- 
pansion is applied directly to the action in ([□]). The expansion should be valid for 
small density of bosons. This is very similar to the discussion of weakly interacting 
bosons |7[]. However, this mean- field solution is not a complete description of the 
physical states of hard-core bosons on a lattice if higher densities are included. This 
can be seen if a nonlocal order parameter is considered. 

b) nonlocal order parameter. Normal Fluid (NF) and Mott Insulator (MI) 

It is possible to relate an order parameter field to the elements of the BWLs by 
replacing w x>x i u X)X > 0, [§]. The partition function then reads 

Nonloc = /exp { - N £ + 2N log det ( y/c + u) }V [u] . (12) 

xx 1 ^x,x' 

Two different contributions can be distinguished, one is the longitudinal component 
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u x ,x+Ae t and the other the transverse component u x>x+e ^. A mean-field approxima- 
tion u X:X+ Ae t ~ °~\\ an d u x, x +e'^ ~ °~± has to satisfy the saddle point conditions 

p- = n --I = 0, P^ = 2da ± -—(n-I) = (13) 
o<j\\ (7|| aa^ a_L 

with the integrals I = a\\ j\ 0(V 2 — Q / ka! p(K)dn and n = j\ 0(k' 2 — Qp(K)dK. 
n is the density of BWLs, k' = 2daa±n + a\\ and p(n) = J 5[(l/d) J2j=i cos kj — 
n\dk\ . . . dkd/(2ii) d . 6 is the Heaviside step function. The mean-field action reads 

s = dal + \<7\\ - log C - J\ log(K /2 /C)e(«: /2 - C)p(«)d« - iir. (14) 
The saddle point conditions C |X3]) imply a 2 + 2^0"^ = 2n. There are three different 



solutions: (<t\\,<t±) = (0,0) (empty phase), (0,^/n/d) (NF) and (V2n,0) (MI). The 
approximation p(n) ~ (1/2)0(1 — |ac|) gives for the mean-field action the values 
shown in Table 1, where ( = 4d(l — n) 2 na 2 for the NF. Although the SF of a) is also 
a phase with fluctuating BWLs, the NF can be distinguished by the fact that its 
fluctuations have a gap in contrast to the gapless fluctuations of the SF. Moreover, 
the SF has a nonzero winding number in contrast to the zero winding number of the 
NF. The fluctuations of the MI also have a gap which is related to an incompressible 
array of BWLs. The latter is commensurate with the lattice. An additional nonlocal 
interaction leads to many commensurate MI phases ||. 

Discussion: For large values of ( and T = there will be always an empty phase 
(i.e., no bosons). (In the case of nonzero temperatures this phase would be a NF 
due to thermal fluctuations.) The parameter a is the hopping rate of the bosons. 
Starting from the empty phase at zero temperature the density can be increased by 
a reduction of (. Depending on the values of a the NF (small a) or the SF (large a) 
will be reached. An increasing density will eventually destroy the SF due to strong 
collisions between the bosons, and a transition to a NF will take place. Finally, 
even the NF will be destroyed: a state can be formed which is commensurate with 
the lattice, provided a is not too large. This is the MI. As an examples the phase 
diagram is shown for d = 3 in Fig.l. 

Table 1: The mean- field solutions for local and nonlocal order parameters 
phase: Empty Phase Fluid Phase Mott Insulator 

density: n = < n < 1 n = 1 

an: V2n 



a ± : ^n/d 

Snonioc: - log C 3n - \ogu - 2 log^v^a) (NF) 1 - log 2 



cp: 2^1 - C/« 4 

s loc : -logC a 2 -(/a 2 -4 log a (SF) 
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Figl: T = phase diagram for hard-core bosons, a is the hopping rate of the 
bosons and y/( the fugacity of empty sites. Thermal fluctuations (T > 0) would 
destroy the empty phase and change it into a NF. 
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